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An analysis is made of the steady-state temperature distribution in a 
poorly conducting plane film on a highly conducting semi-infinite substrate, 
owing to a time -independent heat input in a cylindrical region of the film 
and substrate. The problem is of interest in connection with the localized 
hardening of anodic oxide films on silicon by electron beam bombardment 
in order to produce oxide diffusion masks for the manufacture of integrated 
circuits. A formal solution is obtained for arbitrary dependence of the heat 
input on radius and depth, and a detailed study is made of a particular 
case in which the heat input is independent of radius across the beam, and 
varies in a realistic manner with depth in the film. Approximate formulas 
are given for the temperature in the film when the radius of the beam is 
large compared to the thickness of the film, and also when the conductivity 
of the film is small compared to the conductivity of the substrate. The ap- 
proximate formulas are compared with the results of calculations based on 
the exact solution. Finally, a crude estimate is made of the time required 
to reach the steady state. 

I. INTRODUCTION AND SUMMARY 

Recently considerable interest has developed in the application of 
electron beam technology to microelectronics. 1 A number of papers 
have been concerned with the heat-flow problems encountered when a 
high-power electron beam interacts with a target. Previous investiga- 
tions have considered electron heating of a uniform semi-infinite target, 2 
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of a target consisting of a highly conductive metal film on a less con- 
ductive substrate, 3 and of a thin film not supported by a substrate. 4 
The heating of a poorly conductive film on a highly conductive substrate 
does not appear to have been treated before, and forms the topic of the 
present investigation. It corresponds to the case of an electron beam 
incident upon an oxidized silicon substrate. 

This analysis may find an application in the fabrication of oxide dif- 
fusion masks for integrated circuits by electron beam bombardment. 5 
The etch rate of an anodic oxide film on silicon in hydrofluoric acid has 
been shown to decrease strongly under electron bombardment, and a 
proposal for producing patterns is to harden some areas of the film and 
then to remove the surrounding oxide with dilute HF. 

It should be noted that electron beam bombardment produces radia- 
tion damage as well as thermal effects. The radiation damage alone would 
increase the etch rate of the film in HF, but in conjunction with high 
temperature it also facilitates ionic rearrangement in the Si0 2 film dur- 
ing irradiation, which leads to a decrease in etch rate. The latter effect 
predominates by far in the case of anodic Si0 2 films, so that the net 
result is a strong decrease in etch rate. 

While the rise in temperature during irradiation is thus not the only 
factor contributing to the "hardening" of the oxide film, it is still the 
major factor, and a knowledge of the temperature distribution during 
irradiation is highly desirable. On the one hand one is interested in 
working at a high temperature in order to increase the rate of oxide 
"hardening"; on the other hand one must stay below the melting point 
of silicon (1415°C), or perhaps even lower in order not to generate ex- 
cessive thermal stresses in the silicon. The edge definition of the hard- 
ened region in the oxide film is also of paramount interest for mask fabri- 
cation. Effects due to radiation damage will not be considered here, 
but it may be noted that radiation damage will be generated exclusively 
in the oxide and not in the silicon at the accelerating voltages of interest 
(less than 10 kv). 

In this paper we consider the mathematical problem of calculating 
the steady-state temperature distribution due to an axially symmetric, 
time-independent heat input throughout a cylindrical volume of the 
film and substrate. The thermal properties of both materials are assumed 
independent of temperature, and radiation from the outer surface of the 
film is neglected. A formal solution of the problem is given in Section 
II for an arbitrary dependence of the heat input on radius and depth; 
but in the subsequent analysis we assume that at any given depth 
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the heat input is independent of radius across the beam and zero out- 
side the beam. We also confine our attention to the temperature dis- 
tribution in the film itself. The only thing we really need to know about 
the temperature in the substrate is that its maximum, which occurs 
on the axis at the film-substrate interface, is not high enough to melt 
the substrate. 

With a fixed distribution of input heat, the normalized temperature 
distribution in the film depends on two dimensionless parameters, namely 
the ratio of beam radius to film thickness and the ratio of film conduc- 
tivity to substrate conductivity. In the physical problem, the beam 
radius may be several times the film thickness, and the conductivity of 
the oxide film is between a tenth and a hundredth of the conductivity 
of the silicon substrate. In Section III an asymptotic approximation is 
given for the temperature distribution when the normalized beam radius 
is large. Section IV contains the solution for a perfectly conducting sub- 
strate, as well as an estimate of the first-order effect of finite but large 
substrate conductivity. 

In order to calculate the temperature distribution numerically, it is 
necessary to assume a definite dependence of heat input on distance 
into the film ("depth-dose function"). In Section V we assume a depth- 
dose function which approximates the form determined empirically by 
Griin 6 and also employed by Wells. 7 The parameters are adjusted so 
that the power input is maximum at a depth equal to 40 per cent of the 
film thickness and zero at the bottom of the film, since, in general, 
one wishes to avoid direct heating of the substrate by the electron 
beam. Contour plots of normalized temperature have been calculated 
from the formulas of Section II for selected beam diameters and con- 
ductivity ratios. In addition, the exact temperature distributions along 
the axis and at the top and bottom of the oxide film are compared with 
the approximate formulas of Section III. 

As a typical numerical result, we find that for an SiOo film of thickness 
0.5 micron, bombarded by a 5 kv electron beam of diameter 20 microns 
with a current of 028 ixa and a uniform power density of 10 6 watts/cm 2 , 
the steady-state temperature rise on the axis is about 1800°C at the 
surface of the film, and about 800°C at the surface of the silicon sub- 
strate. 

In Section VI a crude estimate is made of the time required to reach 
the steady-state temperature after the electron beam is instantaneously 
switched on. It appears that in an example such as the preceding, the 
transient time would be of the order of a few tenths of a microsecond. 
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II. STEADY-STATE TEMPERATURE DISTRIBUTION 

The geometry of the problem to be considered is shown in Fig. 1. 
A plane film of thermal conductivity Ki fills the region ^ z ^ c, and 
overlies a semi-infinite substrate of thermal conductivity K 2 which 
fills the region z < 0. We wish to find the steady-state temperature rise 
T{r,z) under the influence of an axially symmetric, distributed heat 
source of strength Q(r,z). 

The temperature rise satisfies Poisson's equation, 



d^T 1 BT , d 2 T _ 
to* r to dz 2 

and the boundary conditions are 



-Q/Kx, 

-Q/K 2 , 



< z < 

z < 0, 



(1) 



T t {r,c) = 0, 
T(r,0 + ) = r(r,<r), 
K 1 7' 2 (r,0 + ) = KiT.irST), 



(2) 



T(r,z)->0 



as 



2 i a 
r + z 



where T z denotes dT/dz. The first of the boundary conditions asserts 
that there is no heat flow across the upper boundary of the film. The 
method of solution which we are going to use would also allow for a 
linearized radiation condition at the surface, i.e., a linear relation be- 
tween T(r,c) and T z (r,c), if one knew the appropriate coefficients. The 
second and third conditions insure the continuity of temperature and 
heat flow across the interface between film and substrate, and the 
fourth condition says that the temperature rise tends to zero at great 
distances from the source. 



V///MW. .;.-. -<r V/M/MM/M 




Fig. 1 — Cross section of plane film on semi-infinite substrate. 
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It will be convenient henceforth to work in terms of normalized, 
dimensionless quantities. In particular, we shall take the film thickness 
as the unit of length, and denote the ratio of film conductivity to sub- 
strate conductivity by e. We also introduce a representative heat 
source strength Qo having the dimensions of power per unit volume. 
Thus, let 

£ = r/c = normalized radius 

f = z/c = normalized depth 

e = KxIKi = conductivity ratio 

q = Q/Qo = normalized heat input 

U = (Ki/cQo)T = normalized temperature rise. 

In terms of these normalized quantities, (1 ) takes the form 

^£4. ld_U,<?U = l-Q 0<f<l, () 

de \ dk df 2 \-eq, f < 0, K ) 

and the boundary conditions (2) become 

U { (Z,1) = 0, 

U(Z,0 + ) = U&CT), 

( 4 ) 
eU^,0 + ) = U { &(T), 

U(U)-+0 as £ 2 + f 2 -*«>. 

In what follows, we shall treat separately the cases of heat input to 
the film and heat input to the substrate. The general case follows by 
superposition. 

2.1 Heat Input to Film 

Assume that q (£,f ) differs from zero only in the film. We denote the 
normalized temperature rise in the film by £/i(£,f), and seek a solution 
of Poisson's equation in the form 

Uiitf) = [ f(w£)Jo(wit)dw, S J £ 1, (5) 

where /(w,f ) is a function to be determined. Similarly, for the normalized 
temperature rise in the substrate, Uz(%,{ ), we seek a solution of Laplace's 
equation in the form 

CTi(W) = f g(w)e w tj (wt)dw, f ^ 0, (6) 

which vanishes as f — * — «> . 
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We shall consider the case in which the normalized heat input can be 
written as 

liter) - *(*)*(r), (7) 

that is, as the product of a function of radius times a function of depth. 
This will probably be justifiable if the increase in beam width with 
depth due to electron scattering is negligible. Furthermore, by taking 
${£) and <p(£ ) equal to 5-functions it is possible to derive the Green's 
function, in terms of which one can express the solution for an arbitrary 
axially symmetric heat input. 

Substituting (5) and (7) into (3) and making use of Bessel's equation, 
we obtain 

[ l/rrKf) - w 2 f(w,t)]J (wl-)dw = -<A(£)<Kf), 

Jo (8) 

for < r < 1. 

From the Hankel inversion formula, 8 it follows that 

f n (w,{) - w 2 f(w,i) = -to* (w)*(f), (9) 

where $ (w) is the Hankel transform of ^(£), denned by 

*M = f #(!)Jo(«*)de. (10) 

Extensive tables 9 of Hankel transforms are available; we note in 
particular the following pairs. 

(i) Uniform beam of normalized radius a: 

(1, 0^i;<«, 

m= \0, S>«, (l la) 

$(w) = (a/w)Ji(aw). 

(ii) Gaussian beam: 

*(t) = exp (-r% 2 ), 

$(w) = (a 2 /2) exp (-aV/4). 
(Hi) Infinitesimally thin, hollow beam of radius £o : 

<K£) = *(€ - 6>), 
$(w) = ZoJo(w£ ). 



(lib) 



(lie) 
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To satisfy the boundary conditions (4), we must have 

/rki) = o, 

f(wfl) = g(w), (12) 

e/ f (tt>,0) = wg(w), 

from which, eliminating g(w), 

fr(w,l) = 0, efs(w,0) = wf(w,0). (13) 

The solution of the two-point boundary value problem for f(w,£) by 
standard methods leads to 



14) 



/Kr) = *<„) Rnh ^r + £ cosh w r f {) cosh w , (1 _ )th 

|_ cosh w -f- e sinn w Jo 

- J <p{ri) sinh io(f - ri)dt\ . 

From the second of (12), 

gM m ^ w) . T <,(„) cosh io(l - ,)d,. (15) 

cosh w + £ sinh w Jo 

In principle, the normalized temperature rise is completely given by 
(5), (6), (10), (14), and (15), provided that the heat input to the film 
can be represented as a product ^(£)¥>(f), and that there is no heat 
input to the substrate. A complete numerical solution for arbitrary 
\p and (p would, however, involve the evaluation of five integrals, each 
of which depends on one or more parameters. In practice, one would 
try to approximate the source function in such a way that at least some 
of the integrations could be done analytically. An example is discussed 
in the following sections. 

2.2 Heat Input to Substrate 

Again we take the heat input in the product form (7), but now assume 
that tp(£ ) differs from zero only in the substrate, f < 0. For the norma- 
lized temperature rise we assume 

Ui(te) = I IW coshw(l - {)J (w£)dw, £ {* £ 1, (16) 
U*(M) = f m(w£)Jo(wt)dw, r ^ 0, (17) 

Ja 
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in the film and substrate, respectively, where l(w) and ra(w,f ) are func- 
tions to be determined. The expression for £/i(£,f) already satisfies 
Laplace's equation and the first of the boundary conditions (4). 

As before, it is easy to show that m(w£) must satisfy the differential 
equation 

mn(w£) - w 2 m(w,{) = -ew$(w)<p(t), (18) 

and the boundary conditions 

l(w) cosh w = m(w,0), 

—el(w) sinh w = m t (w,0), (19) 

m(w£)^>0 as f — > — °°- 

By standard methods we find 

f° 
m(w,S) = D(w)e wi + \e+(w) exp (-w | r, - 1 1) <p(v)dv, . , 

*-* (20) 

for f ^ 0. 

It is clear that ra(w,f ) satisfies the last of the boundary conditions (19) 
if w > and ^(f) vanishes for all sufficiently large negative f. In prac- 
tical cases, it will certainly be justifiable to set the heat input identically 
equal to zero below some finite depth. We shall not take space to in- 
vestigate the mathematical question of how slowly <p(£) could approach 
zero, and still have ra(w,f) also approach zero, as f — > —to. 

From the first two boundary conditions (19), it is straightforward 
to calculate 

l( w ) = *tM , [° exp (wvMv)dv, (21) 

cosh w + € sinn w j-k 

DM = W*) ht W 7 ' ^ " 1 (' exp (wvMv)dv (22) 
|_cosh w -f e sinh wj J-» 

The first of these, substituted into (16), gives the normalized tempera- 
ture rise in the film; and the second, together with (17) and (20), gives 
the temperature rise in the substrate. 

III. APPROXIMATIONS FOR A UNIFORM BEAM OF LARGE RADIUS 

We shall consider henceforth only the case in which the heat input is 
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radially uniform out to the normalized radius £ = a, and zero for £ > a. 
Also we shall be interested only in the temperature distribution in the 
film itself. The dependence of heat input on depth will, however, still 
be taken as arbitrary. 

In the case where heat is applied to the film by means of a radially 
uniform beam which does not penetrate the substrate, the normalized 
temperature rise in the film is given by (5), (11a), and (14). In practice, 
the SiC>2 film may be only half a micron thick while the beam radius is 
several microns. We accordingly seek an asymptotic expansion of the 
temperature distribution for large a. In the analysis we assume that the 
conductivity ratio e is fixed with € < 1 . Our results will also include the 
physically interesting case c« 1. 

Combining (5), (11a), and (14), we may write the temperature dis- 
tribution in the film in the form 



Jo \_w Jq 



Ui(U; e ) = - <p(v)dv + h(w{; e) 



J\{aw)Jn{Pw)dw, 

(23) 

o £ r s i, 

where 

7 / «. \ a Tsinh wt + e cosh wt f , N . ,, x , 

h(w,r, e) = - i— ^ — ; r-r — L / <pbi) cosh w{\ - i))di\ 

w I cosh w + £ sinh w Jo 

(24) 



J <p(rj)dr) — J <p(rj) sinh w(? — r))dr} . 



It is clear that the function h(w,^; e) may be expanded in a power 
series around w = 0; that is, 

00 

h(w,$ ; e) = £ h (m) (0,r, e)w'"/m\, (25) 

7/1=0 

where the superscripts denote derivatives with respect to w. In particu- 
lar, 



(r - e) jf <p( v )d v - jf (r - ,)*(,)d,l. 



(26) 



h(0,r, e) = a 

Let 

p = Z/a, k = ap, (27) 

so that the boundary of the heat input region is p = 1. We have 10 
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Pip) = I Ji(ocw)J (paw) — 
Jo V) 

r dx 

= I Ji(x)J (px) — 

Jn X 

^ p < 1, 

p > 1, 



-E(p), 

IT 



(28) 



*[>©-('-3K): 



where E and K are complete elliptic integrals. 

Furthermore, assuming that p, f , and e are fixed, the following asymp- 
totic expansions for large a are derived in the Appendix, in terms of the 
derivatives of h{w^;e) atw = 0: 



/ h(w,{; €)Ji(aw)J (paw)dw 

Ja 



MQ,r; e) + y (-l)"T(n + J) 



^o r(i)r(n + l) 

s* j-j 4-1.1- ^ ft (2n+1) (0,f; e) 
•F(n + f, » + J; 1; p ) 5-T5 



^ (-l)" +1 r(n + l) 
„^r(i)r(n + i) P 2 «+ 3 



g p < 1, 



p > 1, 



(29) 



where F(a,b; c; 2) is the hypergeometric function. Since when p is near 
unity, we have 



F (n + #, n + *; 1; p 2 ) = 0[(1 - p)" (2n+1) ], 



F (n + I, n + f; 2; i) = 0[(p - l)- (2n+1) ], 



< (1 - p) « 1, 



(30) 



< (p - 1) « 1, 



it follows that the asymptotic expansions (29) are useful only for 
a I 1 — p I » 1, but not in the neighborhood of p = 1. 
Combining (23), (26), (28), and (29), we obtain finally, 



BEAM HEATING OF A THIX FILM G71 



Ui(pa,t; e) ~ caP(p) I vi^d-q + (f - e 2 ) f <p(r,)d v 
- [ (f - vh(v)d V + 0[e/a(l - p)], 



0^ P <1; (31) 
Uiipaf, e) ~ rcxP(p) f <p(r,)d V + 0[e/a(p - Dp 2 ], 

P > I, 

where P(p) is defined by (28). That the remainder terms are 0(f) 
when £ « 1 may be seen from the relations 

h(w£;c) = h(-w,{; -e), 

i(2n+l) / ft v. „\ 1(211+1) / A j. x V -" J " / 

i.e., the odd derivatives of h with respect to w at w> = are odd functions 
of e. Note, however, that setting e = in the asymptotic solution (31) 
does not give the exact solution for a perfectly conducting substrate, 
inasmuch as there are exponentially small terms in a which never appear 
in the asymptotic solution. The exact solution for e = is given in 
Section IV. 

When the product Ea is sufficiently large, the leading terms in the 
asymptotic solution (31) are proportional to P(p); that is, they are 
functions of p (= £/a) only, and are independent of the depth f in the 
film. The function P(p) is plotted in Fig. 2. It is continuous, with a 
logarithmically infinite slope, at p = 1. Numerical comparisons between 
the exact solution (23) and the approximate solution (31) are made in 
Section V. 

We now look briefly at the case of heat input to the substrate by a 
radially uniform beam of normalized radius a and depth dependence 
(pit), for f ^ 0. The normalized temperature rise in the film is, from 
(11a), (16), (21), and (27), 



tt i ,. \ /"" coshw(l — f) r f° , . , ., 

Ui(pa,t;£) = eu I — r- - ^-r- — - / exp {wrj) <p{r))dt) 

Jo w(coshw -\- Esmhw) \_J-co 

•Ji(aw)J (paw)dw, ^ f ^ 1. 



- (33) 



When p, f, and e are fixed, and both a»l and a | 1 — p | 2> 1, an analy- 
sis entirely similar to the preceding gives 
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Fig. 2 — The function P(p) 



-r 



J\{x)J {px) 



dx. 



(34) 



U\(pa,$\ e) ~ eaP(p) / <p(v)dy + e I (q — E)<p(v)dij 

J- oo J- 00 

+ 0[c/a(l - p)], ^ p < 1; 
C7i(pa,f; e) ~ caP(p) / <p(i)dn + 0[e/a(p - Dp 2 ], 

J — co 

P> 1. 



IV. APPROXIMATIONS FOR A UNIFORM BEAM WITH LARGE SUBSTRATE 
CONDUCTIVITY 

We now assume that e <<C 1 ; that is, the conductivity of the substrate 
is large compared to the conductivity of the film. (For an Si0 2 film on 
silicon, e is between 0.1 and 0.01.) Again we consider a radially uniform 
beam, with an arbitrary depth-dose function <piX ). No restrictions are 
placed on the normalized beam radius a. 

For heat input to the film only, the normalized temperature rise in 
the film is given by (5), (11a), and (14). Referring to (14), we expand 
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the function f(iu,£ ) in powers of e. After a little algebra, we find that the 
normalized temperature rise in the film can be written as 

UiiKie) = Eff'lVCfer), (35) 

n=0 ' 



where 



with 



Ux l0 \tt) = a f F(w,t)Jx(aw)J ^w)dw, (36) 



F(iv£) = r— sinh w{ / <p(rj) cosh w(l — n)dri 

w cosh w\_ J; 

+ cosh w(\ — f) J <p( v ) sinh wrj dr\ , 



(37) 



and for n ^ 1 , 



tfi (B) (&0 = (-l) n_1 «/ o fjT ^(,)oodi»(l-,)dfl 



cosh w(l - f) tanh w , , w ,. . dw 
cosh 2 w w 



(38) 



The quantity t/i (0) (£,17) is the normalized temperature rise in the film 
when the substrate is perfectly conducting. In this case, the temperature 
rise at the bottom of the film is zero, and in fact it is obvious from (37 ) 
that F(w,0) = 0. The quantity ?7i (B) (£,$■) represents the nth order 
correction if e is small but not zero. 

We may evaluate f/i (0) (£,f ) by contour integration. Let S„ denote the 
semicircle of radius nir in the upper half-plane (n = 1, 2, • • •), with 
diameter along the real axis indented at the origin. From (37) it follows 
that F(w,t;) is uniformly bounded on S n , and has simple poles within 
S„ at w = (ra — |)W (m — 1, • • • , n). The choice of integrand depends 
on whether ^ £ ;£ a, or £ g£ a. For ^ f 5j a we consider 

f F(w,i;)MZu))Hi m (otw)dw, (39) 

and for £ ^ a we consider 

f F(w,t)Ji(aw)H (1) ^w)dw, (40) 

where H Q W and Hi W are Hankel functions. The integrals are evaluated 
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by the calculus of residues and then the limit n — > °o is taken. Since the 
procedure is a standard one, we omit the details and merely state the 
results. 

We obtain, for ^ £ ^ a, 



Ui w (U) ~ [f / <p(v)dv + I w(v)dv 



= _ 2a £ sin Km ~ \U] h[(m _ hM]Ki[{m _ i)va] (41) 

• / ip(n) sin [(ro — |) TJjlf/T?, 
•'o 

and, for £ ^ a, 

( 0), . .» 2a A sin [(m - f )*£] 

t/i KM) = — L. / _ n 

tt m =i (m 2 ) (42) 

•/i[(m - i)7ra]K [(w - §)*•£] / <^(r?) sin [(m - %)ni\dn. 

Here h , h , K , and Ki are modified Bessel functions. The continuity 
of dUi 0) /d£ at £ = a is readily verified, while that of £/i <0) at £ = a 
follows from the identities 

Io(x)Ki(z) + 7i (*)*#(*) - 1/*. ( 4 3) 

and 

E C0S , [(w "* M =^(l-g), 0g^2. (44) 

£=i (m — j) 2 2 

We remark that (41) and (42) could have been derived by separation 
of variables in (3). For | a - £ | » 1 the right-hand sides of (41) and 
(42) are exponentially small. It follows that (41) and (42) are con- 
sistent with the asymptotic expansions (31) and (32) for a. » 1, when 
c = and asymptotically small terms are neglected. 

It does not appear possible to evaluate the first-order correction 
Ui ll) (M)> as given by (38), using contour integration, because the 
integrand has the wrong parity in w. We can, however, obtain a bound 
on the value of U\ l) (0,0), at the "hot spot" of the film-substrate inter- 
face, where the zero-order solution Ui (0) (0,0) vanishes. 

From (38), setting n = 1 and changing the order of integration, 
which is justified since the double integral is absolutely convergent, 



C/i (1) (0,0) -a f\(v)\j 
•'o •'o 



cosh w(l — 7j ) Ji(aw) 7 

; aw 

cosh w w 



d v . (45) 
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When 77 = 0, the inner integral is equal 1 " to 1. When 77 > 0, we trans- 
form the inner integral by substituting the integral representation 13 

rx/2 

Ji(aw) = (2aw/ir) I cos (aw cos 0) sin 2 dd, (46) 

•Jo 

and again invoking the absolute convergence of the double integral to 
change the order of integration. This leads to 



f. 



cosh w(\ — 77) Ji(aw) j 

-. aw 

cosh w w 



2a r'TrcoshwU -77) , flU 1 . 2 ,. 

= — / / -. cos yaw cos 6)clw sin dd 

it J o \_ J o cosh w J 

. 7T7? f' 2 cosh (^Tra cos 0) sin 2 ^0 ,._ N 

= a Sin — / . , , (X tt— ; rryi — r (47) 

2 Jo sinn- \\ica cos 0) + sim (57777) 

. 7T7j /" x/2 cosh (\ira cos 0) sin dd 
2 Jq sinh 2 (^7ra cos 0) -f sin 2 (57717) 

2 _i fsinh (f7ra)~| 
= - tan . ' — r < 1, 

where the third line follows from a table of Fourier transforms. 14 Hence, 
finally, from (45) and (47), 

c/! (1) (0,0) 

2 (' , v . irr/r/** 72 cosh (^Tra COS 0) sin 2 d dd ~|, 

< a / p(»j)dt7. (48) 

We see from (31) that asymptotically, for a ~» 1, the upper bound in 
(48) is attained, since P(l) = 1. 

Now suppose that heat is put only into the substrate, so that the 
normalized temperature rise in the film is given by (33). Expanding in 
powers of £ leads to 

Udtf;e) = Zf n r 1 ,n> a,f), (49) 

n = l 

where 

( - 1 ) " _1 cosh (1 - f ) w tanh" -1 w 



•'0 



w cosh w 



X / exp (wn)<p(rf)drj \Ji(aw)J (£w)dw. 



(50) 



(51) 
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In particular we have, on changing the order of integration, 

J—ao _Jq W J 

after substituting the known value 15 of the inner integral. 

V. NUMERICAL RESULTS 

In this section we give the results of some numerical computations 
using the exact formulas of Section II, and some comparisons with the 
approximations of Sections III and IV. We asume that the electron 
beam voltage is such that the electrons penetrate to the bottom of the 
oxide film (about 5 kv for an 0.5 m film 7,16 ), but do not enter the sub- 
strate. For the depth-dose function we take 

p(£) = sin #, = 5ir/6, £ t £ 1, (52) 

which is plotted in Fig. 3. The assumed depth-dose function vanishes at 



< 0.4 



0.2 



0.5 

HEAT INPUT, 0(f) 



Fig. 3 — Heat input function <p(f) = sin 5tt{-/6. 
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the bottom of the film, has a maximum at a depth equal to 40 per cent 
of the film thickness, and in general corresponds very closely to the 
empirical function used by Grim 6 and Wells. 7 The position of the maxi- 
mum could be varied, of course, by changing the parameter /3. 

Substituting tp(£ ) into (5), (11a), and (14), we find after some al- 
gebra that the normalized temperature rise in the film is given by 



UdU) = aV(i)Bmft + a0W(te), 



where 



and 



vis) = f 

•In 



Ji(atw)J (£w) 
o w 2 + /3 2 



dw 



(53) 



(54) 



W(tf\ [" Ji(aw)Jo(£w) 

U ' n Jo w{w 2 + $ 2 ) 

Ye cosh(l — f )w — (sinh f w + e cosh $w) cos /3 



(55) 



cosh w + £ sinh w 



dw. 



The integral on the right side of (54) can be expressed in terms of 
modified Bessel functions. We have 1 



.(, 



wJi(aw)Ji(tw) 
o (w 2 + /3 2 ) 



dw = 



\li(pt)Ki(fia), 
[li(0a)KiW, 



t g a, 
t ^ a. 



(50) 



Integrating both sides with respect to t from a to £ and using the rela- 
tionship 18 



/ 



dw = 



o w 2 + /3 2 



(57) 



we obtain 



Ftt) = 



j8 Jot 

h($a)K a (M) 



]■ 







^ £ ^ a, 



(58) 



The integral for TF(f,f), on the other hand, has to be evaluated 
numerically. The integrand is oscillatory, and falls off exponentially for 
large w if < f < 1. If f = or f = 1, it falls off like 1/w* if £ ^ 0, 
and like 1/w 2 if £ = 0. The numerical integration was done by Simp- 
son's rule on an IBM 7094 computer. Combined analytic and empirical 
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investigations of the accuracy were made in order to guarantee that the 
relative error in any value of U is less (in most cases, much less) than 
one per cent. 

Four different sets of parameters were chosen; namely, a = 2, 10, and 
20 with c = 1/40, and a = 2 withe = 1/4. The normalized temperature 
rises at the surface of the film and at the film-substrate interface are 
plotted against normalized radius in Fig. 4. Note the differences in 
scale; in each case the edge of the beam, £ = a, is at the center of the 
plot. The temperature distribution along the vertical axis is shown for 
the same four cases in Fig. 5. 

It is seen that the temperature distribution at the surface of the film 
becomes more flat-topped, and the fall-off at the edge of the beam 
becomes relatively (although not absolutely) mote abrupt as a increases 
in the first three cases of Fig. 4. Also note that the temperature levels are 
somewhat higher and the temperature variation through the film is less 
in Fig. 4(d) than in Fig. 4(a), since for the same value of a the relative 
conductivity of the substrate is only xV as large in Fig. 4(d) as in Fig. 
4(a). 

The dashed curves in Figs. 4 and 5 correspond to the approximate 
formulas (31 ) for large a. If *>(£) is given by (52), these approximations 
read: 

,r / -, r>/ \ 1 — cos 6 
Ui(pa£; e) ~ eaP(p) ? 

sin j8f 
L~T~ "J' (;i!n 

£ p < 1, 



+1 



- £ cos/3 - £'(l - cos/3) , 



,, / \ tw \ 1 — cos 8 . , 

Ui(pa,r, e) ~ caP(p) * , P > 1, 

where p = £/a. As expected, the approximations are discontinuous at 
the edge of the beam, p = 1 ; and they are not much good when a = 2 
(worse for the larger value of e). They are remarkably good, however, 
for a = 10 and a = 20; the dashed curves essentially coincide with the 
solid ones except on the surface of the film in the immediate neighbor- 
hood of the beam edge. 

Contour plots for the temperature distribution in the film are given 
in Fig. 6 for a = 2 and a = 10 with £ = 1/40, and for a = 2 with 
e = 1/4. Contour plots were not made for a = 20, because the numerical 
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Fig. 4 — Normalized temperature rise at upper and lower surfaces of film: 
(a) a = 2,e = 1/40; (b)a = 10, c = 1/40; (c) a = 20, e = 1/40; (d) a = 2, e = 1/4. 

integration is slow for large a (the integrands oscillate more rapidly); 
but it is clear that for large a the approximate formulas (59) would 
yield accurate contours, except very close to the beam edge. 

We may also compare the bound on the first-order correction term 
for small e, as given in Section IV, with the exact results. At the center 
of the film-substrate interface, (48) and (52) give 

ftV"(0,0) < e«(l - cos/3)//3 = 0.713**, (60) 

which leads to the following comparison with the exact solution Ui(0,0). 



680 THE BELL SYSTEM TECHNICAL JOURNAL, MAY-JUNE 1966 



»*> 0.2 



N 1.0 





(a) 


u(o,n/ 


i 
i 
i 

i 
i 


EXACT 

APPROXIMATE 


/ / 
/ / 


i 
i 
/ 




a = 2 
* = 40- 




/ / 

r / 






X' 










/' 









1.0 



0.6 



0.2 





(b) 






u(o,n/ 
















a = io 














./ 


/ 






> 


/ 











(c) 










u 


My 




O = 20 

-A 























0.4 



0.2 



4 0.6 0.8 0.2 0.4 

NORMALIZED TEMPERATURE RISE, U 



Fig. 5 — Normalized temperature rise along axis: (a) a = 2, c = 1/40; (b) 
a = 10, £ = 1/40; (c) a = 20, £ = 1/40; (d) a = 2, £ = 1/4. 
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In order to relate the dimensionless temperature rise U\ to the physical 
temperature rise 7\ , it is convenient to introduce the power density 
(i.e., per unit area) in the incident beam. For a uniform beam of radius 
ac with depth -dose function <p(z/c), the dimensional factor Q , which 
normalizes the heat input per unit volume (Section II), is related to 
the incident power density Po by 



0.1 0.2 0.3 0.4 0.5 0.1 0.2 0.3 0.4 0.5 0.€ 

1.0 
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ir(ac) 2 Po = 2irQo I / <p(z/c)r dr dz, 
Q = Pa/cj fp(f)df . 
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(61) 
(62) 



Hence, the actual temperature rise at the point (r,z) of the film is given 
by 

Tl (r,z) = (cWKi)tfi(r/c, z/c) 

(63) 
= (1.403cP /^i)?7i(r/c,z/c), 
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Fig. 6 — Isothermal contours, (a) a = 2, e = 1/40; (b) a = 10, e = 1/40; (c) 
a = 2, e = 1/4. 
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where the numerical coefficient corresponds to the depth-dose function 
(52). Consistent units for (63) are: 

Ti = temperature rise in °C 
r,z = coordinates in cm 

c = film thickness in cm 
Po = incident power density in watts/cm" 
Ki = film conductivity in watt/cm °K 
0.239 Ki = film conductivity in cal/sec cm °K. 

Wc shall now look at a numerical example. It turns out that the 
constant-conductivities model which has been analyzed in the present 
paper is not a very good approximation to the real problem of a silicon 
dioxide film on a silicon substrate, because the thermal conductivities 
of both materials depend strongly on temperature. In fact, the conduc- 
tivity 19 of SiC>2 increases from about 0.015 watt/cm °K at room tempera- 
ture (various values are reported — not for thin films — which differ 
among themselves by as much as 2 : 1 depending on the crystalline orien- 
tation of the sample) to about 0.03 watt/cm °K at 900°C. For Si, on the 
other hand, the conductivity decreases from about 1 watt/cm °K at 
room temperature to about 0.03 watt/cm °K at 900 °C. For purposes of 
calculation we shall more or less arbitrarily assume the values 

Ki = 0.03 watt/cm °K 
Ki = 1.2 watt/cm °K 
c = 1/40 (64) 



\— 5 



c = 0.5 m = 5 X 10 cm 
Po = 10 watts/cm . 

Since these conductivities may be somewhat larger than the actual 
conductivities, the temperatures which we shall compute may be some- 
what lower than the actual tempera tmes. For a 5-kv beam, the assumed 
power density corresponds to a current density of 200 amps/cm 2 . 

Table I gives the total beam current and the maximum temperature 
rise (i.e., on the axis) at the top and bottom of the film, for beams of 
diameter 2 n, 10 n, and 20 n, corresponding to the previous computations 
with a = 2, 10, and 20, and a conductivity ratio of 1/40. It appears, 
therefore, that in each case at least a part of the irradiated spot would 
be raised to the temperature at which the oxide hardens (900°C), but 
in no case would the substrate melt (1415°C). 

It is probably worth repeating that the physical problem of interest 
is nonlinear, because of the dependence of conductivity on temperature. 
Bounds on the solution may be obtained from linear models, by using 



BEAM HEATING OF A THIN FILM 

Table I 



083 



Diameter 


Current 


Ui(0,l) 


l/i(0,0) 


ri(o,c) 


Ti(0,0) 


2m 
10^ 
20/i 


6.28/xa 
157/ia 


0.3857 
0.5799 
0.7589 


0.0314 
0.1771 
0.355G 


902° 
1356° 
1775° 


73° 

414° 
832° 



the theorem that with a fixed heat input the steady -state temperature 
is not increased anywhere (usually, it is decreased everywhere) if the 
conductivity is increased anywhere, and vice versa. However, only a 
full-dress numerical treatment of the nonlinear partial differential 
equation, assuming that one knew the temperature dependence of the 
conductivity, would be likely to yield really accurate results. 



VI. TRANSIENT EFFECTS 

It is of interest to know how long it will take to reach the steady 
state if the electron beam is suddenly switched onto the film, since this 
gives an idea of how rapidly the beam may be scanned in laying out a 
mask. There have been some published analyses 7 - 21 of transient heating 
effects in electron beam machining, but we shall content ourselves with 
a crude estimate of the time scale in the present problem. 

Consider the case of a film on a perfectly conducting substrate, with 
the film initially at zero temperature, and with a time-independent heat 
input starting at t = 0. Then the instantaneous temperature distribution 
satisfies the heat flow equation 



tfT , 1 dT 

dr' 1 r dr 



d'T 
dz* 



Q_,Ci9T 
Kx K x dt 



(65; 



where K x is the thermal conductivity, C the heat capacity, and 8 the 
density. The total temperature T(r,z,t) may be written as the sum of a 
steady -state part and a transient part, 



T(r,z,t) = T,{r,z) + G(r,z,t), 



(00) 



where Ti(r,z) satisfies Poisson's equation (cf. Section II) and Q(r f z,t] 
satisfies the homogeneous equation 



ye , l ae , a 2 e cs ae _ Q 

dr- r dr dz n - Ki dt 



(67) 



and vanishes as / — * <x> . 

A sufficiently general solution of (07) may be written in the form 
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e(rvM) = E r A m (w) exp [-(Ki/C6)[w* + (to + $)V/c 2 ]t} 

m=0 •'O 

. . r / v . (f» + £)*« , 

X Jo(wr) sin 2 tfw. 



(68) 



The initial condition requires that the total temperature vanish at 
t = 0; that is, 

T,{r,z) 4- 6(r,z,0) = 0. (69) 

Hence, from a knowledge of the steady-state temperature one can in 
principle use the properties of Fourier series and Fourier-Bessel integrals 
to determine the functions A m (w) and the transient solution Q(r,z,t). 
We would like to know how fast Q(r,z,t) approaches zero with in- 
creasing time. It is clear that (68) cannot be characterized by any single 
exponential decay; but we observe that the most slowly decaying ex- 
ponential is 

exp [- (iLV/4c 2 C5)f], 

and it is therefore reasonable to define a crude "transient time" as 

r = 4c 2 C5/x 2 X 1 . (70) 

We assume the following numerical values for the SiC>2 film: 



c =5X10 cm 
C = 1 watt sec/gm °K 

5 = 2.2 gm/cm 3 
K x = 0.03 watt/cm °K. 



(71) 



Then 

v- 7 



r = 3.0 X 10 ' sec, (72) 

so the transient time is a fraction of a microsecond. 
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APPENDIX 



Asymptotic Expansion of a Class of Integrals 
In this appendix, the asymptotic expansion of 

J(p,a) = I h(w)Ji(aw)J (paw)dw (73) 

•'o 

is derived, where 

00 

hiw) = ^h (m) (0)w m /m\, (74) 

ro=0 

for a >>> 1 and a | 1 — p | ^> 1, that is, for p not in the neighborhood of 
1. It is clear that the asymptotic expansions will break down in the 
neighborhood of p = 1, since then the integrand will contain a term 
which is not rapidly oscillating. 

We start from a result given by Tranter. 22 Namely, if we have the 
expansion 

i»00 00 

/ exp (— yw) 5(a,w)dw = 2 A n (a)y n , (75) 

•>0 71=0 

then, formally, 

[ h(w)5(a,w)dw~J2 (-l) n A n (a)h (n) (0). (76) 

•'0 n = 

In the case at hand, 

5(a,w) = Ji(<xw)Jo(paw). (77) 

Assume first that p is fixed and ^ p < 1. Then if (a 2 + y 2 )* > pa + y, 

1 ^ (-i)V"r(2m + 2) 
'o 2 2m + 1 [V(m + l)] 2 

/ 2 \m+l / 2 \ ( 78 ) 

Using standard transformations, 2 



we have" 



rl 
exp (yw) J 1 (aiv)Jo(paw)dw = - V 
d tn=( 
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.( 



exp (—710) Ji(aw)Jo(paw)dw 

1 A (-l)"'p 2m r(2m + 2) /«Y +1 
a^o 2 2 "'+"[r(m + l)] 2 \7V 

XF(ffl+ l,m + f;2; - ^J 

1 , 7 y> (-l)V m r(2m + 2)r(-j) 
a a 2 & 2 2 '"+ I [r(w + l)] 3 r(| - m) 

XF(m+f,m + i;f; -^ (79) 

= 1 _2t_ A P 2 "T(m + j)r(m + 1) 
a Tra 2 ^ [r(m + l)] 2 

X^H,m + i;!;-J 

I + _J_ T (-D" +1 7 2n r(n + i) 
a VW& a 2 T(w + l) 

X F(n + t,n + i;l;p 2 ), 

where in the last step we have expanded the hypergeometric function 
in a power series and interchanged the order of summation. Comparing 
(79) with (75), we see that 

fl/a if n = 0, 
A 2n ( a ) = J (80) 

[0 if «> 0, 

A I \ ( _ D" +I r(n + i) „/ , 3 , i t 2x 

iWo) = a ^ra)r(»+i) n« + i,« + i;Up). 

It follows from (76) that 
/ h(w)Ji(aw)Jo(paw)dw 

h(0) , A (-l)"r(n+j) , v 

~ — + & r(«r(n + i) (81) 

XF(» + !,n+l;l;p ! )^F 

a 2n+2 

for < P < 1. 
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An entirely similar derivation, the details of which will be omitted, 
leads to the expansion 

(-i) n+1 r(n + f) 



/ h(w)Ji(aw)J (paw)dw ~ 2^ r m r / . n ,„ 
Jo m=o I [3)1 [n -+- L)pr 

forp > 1. 



(82) 
} (0) 

2n+2 ' 
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